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Investigation of many problems of acoustics, hydromechanics, elastodynamics and others 
C\) | is connected with boundary value problem (BVP) for wave equation - the multidimensional 
analogue of the Dalamber's equation, which describes processes of spreading the waves in uniform 
isotropic ambience. So elaboration of the efficient methods of their solving for areas with free 
geometry and varied type of the border conditions is very currently . 

The most efficient method of the study of such problems is the boundary integral equations 
method (BIEM). Its main advance is reduction of dimension of solved equations and increasing 
PLh ' calculation stability. That is very greatly, particularly for unbounded areas. The existing meth- 
ods and programs of splines approximation for free contours and surfaces remove the problem 
of restriction of the form of considered areas by using BIEM. At present time this method is 
broadly used for solving elliptical and parabolic problems, what is connected with success in 
development of BIE theories for elliptical and parabolic equations and systems. 

Solving nonstationary dynamic problems on base of the BIE method requires entering the 
notion of the generalized solution. That is connected with particularity of the fundamental 
solution of hyperbolic equations, which belongs to the class of generalized functions. Besides 
■^J- ! classical notion to differentiability of the decisions for hyperbolic equations sharply narrows 
the class of useful for applications problems. In particular, typical physical processes, being 
accompanied shock waves, are not described by differentiated solutions of such equations. 



Here BIE method is elaborated for determination of nonstationary solutions of wave equations 
in bounden region in R N with boundary Dirihlet or Neumann conditions. The generalized 
solution of BVP are constructed on the base of dynamic analogues of Green and Gauss formulas 
for solutions of the wave equation in space of generalized function. Their regular integral repre- 
sentations and singular BIE are built for N=l,2,3, also at presence of shock waves. 
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1. The generalized solutions of the wave equation. Shock waves. The multidimen- 
$2j ■ sional analogue of the Dalamber's equation is considered: 

1 d 2 v 

n c u = Au--— = G(x,t), xeR N } teR\ (1) 
& ot z 

G is regular function, here and hereinafter A is Laplace operator, u,i = J^-, ii = u, t = % . 

It is well known [1,2] that Eq. JH) is strictly hyperbolic, class of its solutions contains functions 
with breakup of derivatives. The surface of the breakup F in R N+1 = {(x,r = ct)} is the 
characteristic surface of Eq. (fTJ) , which satisfies to the relations: v 2 — ||z/||^ = 0, v T = i>n+i < , 
where u(x, r) — {v\, v^, v T ) is a normal vector to F , \\v\\ N = ^VjVj . In R N this is a wave 
front F t , moving with the constant speed c : 

c=-v t /\\v\\ N , v t = cv T (2) 

On reiterative indexes i,j in product everywhere their is summation from 1 to N . 
There are the Adamar's conditions on the jumps at F t : 

[u (x, t)] Ft = 0, [u + cm u,i] Ft = 0. (3) 
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Here [f(x,t)] F = f + (x,t) — f (x,t) — lim (f(x + en,t)—f(x — en,t)), x G F t , n(x,t) is 

* ' e— >+0 

a unit normal vector to F t , directed aside its spreading: 

Ui gradF . 

ni = Tir = T\ i = h-,N; 4 

\\u\\ N \\gradF\\ 

Last equality may be written only if the equation of wave front can be present as F(x, t) — , 
at condition of existence grad F . 

The class of the similar solutions of the hyperbolic equations is named shock waves. 

From second condition Q follows 

iT + crii u~ = ii + + crii u,l , x G F t . (5) 

If before the front u = , this equality gives the useful correlation on F t : (gradu, n) = —c~ l u. 

We notice, that tangent derivatives to characteristic surface, on the strength of continuity 
u , are also continues. Then 

[u, T It + uftUp = for V 7 G R N+1 : {v, 7) = 0. (6) 

In particular if 7 = 7 J = (— Vj, v T 5{, v T b\, u T 5 3 N ) , where b\ is Kronekker's symbol, then we 
have the condition of the type: 

[-u, T Vj + UjU T } F = [wnj + ratj-]^ = 0, j = 1, N. (7) 

Hereinafter we shall consider the functions u(x, t) , which are continues together with derivati- 
ves of first and second order almost everywhere with the exclusion of finite or counting number 
of wave fronts, on which the conditions on jumps are satisfied (J2j) . We shall name such solutions 
classical. Let's show that they are generalized solutions of the Eq. JTJ). 

For this let's consider the Eq.JTJ on the space of generalized functions D'(R N+1 ) = 

= r)|, which are determined on the space of infinitely differentiable finite functions 

D(R N+1 ) = {ip(x, t)} [2]. The value / on ip , as it's accepted , is denoted as (/, if) . For regular 
function / , corresponding to local integrable / , (/, ip) = f f(x, r)ip(x, r)dV(x)dr . Here 

R N+1 

everywhere dV(x) = dx\...dxN ■ 

Definition. Function / G D'(R N ) is identified as generalized solution of Eq. JTJ if for 
any up G D(R N+1 ) (□„/, y) = (/, D c (p) = {G,cp). 

Lemma 1.1. // u(x,t) is the classical solution of Eq. (QP, then u(x,t) is its generalized 
solution . 

P r o o f. If u(x, t) has a finite breakup on F then [2] u,j = u,j +[u]pVjSp(x,T), where first 
summand on the right is a classical derivative over Xj, \\u\\ = 1, 5f is simple layer on F : 



([u) F UjS F (x,r),<p(-K,r)) = / [u(x,r)) F fj(x,T)<p(x,r)dF(x,r), Vcp G D(R 



F 



Here integral on F surfaced. On the strength of continuity u outside of wavefront 
[u] F = lim (u(x + en, t) — u(x — £7i,t)) = u + (x,t) — u~{x,t) = [u] F . 

So, with regard for (pjj), we get: 
u,j = u,j +[u] F u j 5 F (x,T), u,jj = u,jj + [u,j] F \\v\\ N nj5 F + d,- {Hjv Mf^f] ■ 



On virtue of Eq.Q 

U, T = U, T +[u} Ft U T 5 F = C~ l U, t - \v\ N \u\ F b F , U, tt = C~ 2 U, tt ~C~ X [u, t ] Ft \\u\\ N S F ~{ \\u\\ N [it] ^ Sp } ^ . 

With regard for these equalities and Adamar's conditions (|T]) 

D c u = D c u + {cT 1 [u,t] Ft + [njU,j) Ft } \\u\\ N 5 F (x,r) + 
+c~ 1 d t {\\v\\ N [u} Ft 5 F (x,T)} + dj {\\u\\ N [u] F n j 5 F {x,r)} = G(x,t), 

since all densities of simple and double layers on F t are equal to zero. Really, the second 
summand is a zero on the strength of the second condition Q on front. But two others are null 
on the strength of the first one, since their action on D(R N+l ) is defined as 

(c~ l d t {\\u\\ N [u] Ft 5 F (x,r)} + dj {\\u\\ N [u] Ft njS F (x,T)} ,<p(x,t)) = 

= ~ J v T [u) F - yv J dF(x, r) = 0. 

F 

The lemma has been proved. 

Remark 1. From this lemma follow that condition on fronts of the shock waves easy to 
get, considering classical solutions of the hyperbolic equations as generalized one. It is enough 
to equate to zero the density, corresponding to independent singular generalized functions - 
analogues of simple, double and others layers, appearing under generalized differentiation of 
the solutions. The determination of such conditions on base of the classical methods is more 
labour-consuming procedure. 

Remark 2. Eq.(JH) allows to consider the generalized solutions with derivatives breakup 
also on moving surface F(x, t) — , which velocity of the motion can depend on point of the 
front, then v(x,t) = — F, t / \\grad F\\ . On they the Adamar's conditions @are fulfilled with 
change c to v(x,t) . Such solutions can be generated by the right part of equation if the support 
of function G(x, t) enlarges in the course of time in R N . 

2. Statement of nonstationary BVP. Uniqueness of the solution. Let the solution 
of Eq.(|TJ is determined in S~ C R N , bounded by Lyapunov's surface S ([2], p. 409), t > . 
Initial condition: by t = 

u(x, 0) = u (x) for x G S~ + S, u, t (x, 0) = u (x) for x E S~ . (8) 

We'll consider two BVP problems, consequently Dirihlet or Neumann boundary conditions: 

u = u s (x, t) for x G S, t > (first BVP); (9) 



— = p(x,t) for x G S, t > (second BVP). (10) 
on 

Here we denote as n = (ni, ...,un) a unit vector of external normal to S , D~ = S~ x R + , 
R + = [0, oo) . 

It is supposed that initial and boundary functions Uq(x) and us(x,t) are continues, uq(x) 
and p(x, t) piecewise continues. For the first BVP the boundary and initial conditions are 
coordinated: 

Uq(x) = us{x, 0) for x G S. (11) 
On wave fronts, if they appear, Adamar's conditions (JBJ) are fulfilled. 
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Notice that shock waves always appear if there is not coordination condition between initial 
velocities and velocities at boundary: 



Uq(x) = its(x, 0) for x G S, 



(12) 



That is typically for physical problems. In this case at initial moment of time at the boundary .D 
the shock wave front is formed, which spreads with the speed c in R N . For building continuously 
differentiated solutions this condition is necessary. Here we don't suppose it. 

It's supposed that initial conditions have been given and the one of the boundary conditions 
is known accordingly considered problem. 



N 



N 



We enter the functions E = 0, 5 \u 2 + M >? > ^ = ^> ^ M >? — S u 



3=1 



3=1 



Lemma 2.1. If u(x,t) is a classical solution of Eq. J2JJ then 



[E] 



F t 



— C 



u 



du 
dn 



[L(x,t)] F = c 2 [u + c 



F t 



du 
dn 



u . 



(13) 



Proof. On the strength of equality [ab] = a + [b] + b [a] , with regard to JTJ and 0, we 
have 



„ du 
cE + u— 
on 



du 



0, 5 (c u 2 + cu,jU,j) + u— 

on 





du 




u + c— 
dn 



+ (u + cu,j Uj) [u)J + 0, 5u,+ [ cu,j +uuj] + 

+0, 5 (cu,j +u~nj) [u,j ] = 0, 5c _1 [u] (u~ + cu~ n,) + 
+0, 5 [u,j] (cu~ +u~nj) = 0, 5cu,j [u,j +c~ 1 nju] + 0, 5c~ x u~\cnjU,j +u] = 

(here n is a unit normal to front F t in R N ). Thence the first formula of the lemma follows. 

Hereinafter since [a 2 ] = (a + + a~) [a] , and on the strength of Q and ©, we get second 
formula (HI 



[L] = 0,5 



N 



u. 



3=1 



0, 5 (u + +u, T ) [u, T ] -0,5 (u,t ) [ U)i ] 



0, 5 (w + +u, T ) [«, T ] +0,5 («,t +u, • ) [ tt, r ] 



0, 5 { (m,+ + (u, T ) } [ u, r ] = c 2 iii +c 



du 
dn 



u 



Remark. If before the front of the wave u = , that, with use of (jHJ), we have [L(x, = 
i.e. in this case function L continues. 

Theorem 2.1. // u(x, t) is classical solution of BVP, then 



(E{x,t) - E(x,0))dV{x) 



t t 
dt / G(x,t)u, t dV(x) + / / us(x,t)p(x,t)dS(x)dt. 




o d- 



o s 



Proof. Multiplying Eq. ((TJ) on u, T in the field of differentiability, after simple transforma- 
tions we get: 

E, T -(u, T u,j),j= -u, T G. (14) 



4 



But now we integrate Eq. ljTlj) on D , with regard to dividing the area of integration D on 
the parts, bounded inside by wave fronts F k . 

Let us consider the left part of this equation as divergency of some vector in R N+1 , which is 
continues and differentiable in region between fronts. Then, using Ostrogradskiy-Gauss theorem 
in R N+1 , we get 



E, T dV(x)dr — j (u, T u,j),j dV(x)dr + J u, T G(x,r)dV(x)dT 

D- D~ D~ 



u, T G(x,T)dV(x)dr + / (E(x, r) — E(x, 0))dV(x) — / / (u, T u,j %) dS(x)dr+ 




TP. ^ 



- U, T U,j Vj] Fk dF k (x, t) = 0, 



here dF k (x,r) is a differential of the area of surface in corresponding point of wave front F k . 
On the strength of formulas ^ and (fTTH) . \Ev T — u, T u,j Vj] F = —c^W^Wj^ [cE + Uj£\ = . So 
the last integral is equal to zero. With regard for indications for boundary function, thence the 
formula of the theorem follows. 

Remark. First condition of formulas (fTTf easy be get, considering Eq. ((H)) in D'(R N+1 ) : 
E, T -(u, T u,j),j= —u, T G + {[E] v T — [u, T u,j)vj}5 F = -u, T G- \\u\\ N {[E) + [u^u^rij} 5 Ft . 
Eq. (USD is fulfilled in D'{R N+1 ) only if [E + c^u^]^ = 

Theorem 2.2. // classical solution of the first (second) BVP exists, then it is single. 

Proof. On the strength of linearity of the problem, it is enough to prove uniqueness of null 
solution. For it G = , initial and corresponding boundary conditions are null. Then, as easy 

see, from theorem 2.1 follows that J E(x,t)dV(x) = 0. Since E is nonnegative, consequently 

s- 

E = =3- u = const. From initial conditions follows u = 0. 

T e o r e m 2.3. If u(x,t) is classical solution of the BVP, then 



t 

L(x,t)dV(x)dt = J uG(x,t)dV(x)dt — J J us(x,t)p(x,t)dS(x)dt 



d- d- os 

-2 



+ 

+c~ 2 J (uii(x, t) — Uoiio(x)) dV(x) 
s- 

Proof. Multiplying Eq.JTJ to u , after simple transformations we have: 

L(x, r) + (uu,j ),j —{uu, T ), T = G. (15) 
Let integrate (jTBT) over D~ with regard for its partition with wave fronts F k . Similarly, as in 
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theorem 2.1, using Ostrogradskiy-Gauss theorem , we get 



LdV(x)dr = J uGdV(x)dr + J ((uu, T ), T —(uu,j),j) dV(x)dr = 

L> D~ D- 

T 

uG dV(x)dr + c _1 J (uu, t —u u ) dV(x) — J dr J uu,jUjdS(x)+ 
d- s- OS 



+ / [uu, T p T -uu, j i'j] Fk dF k (x,T) 



On the strength of Adamar's conditions JH|),the last summand is equal to zero. So, with regard 
for conditions on boundary ((HJ), (fTTIll . we get the formula of the theorem. 

Remark. Theorem 2.1. is the law of energy conservation. Similarly it's convenient to use 
such law with regard for conditions on front of the shock waves for proof of uniqueness of BVP 
solutions in nonstationary problems of mathematical physics . 

3. The dynamic analogue of Green formula in D'(R N+1 ). For building the BVP 
solution we consider it in D'(R N+1 ) . For this we introduce the characteristic function of semy- 
cylinder D~ : H^(x,t) = Hg (x) H(t) , where Hg (x) is characteristic function of set S~ , which 
is equal to 0,5 on 5* , H{t) is Heaviside function equal to 0,5 by t = . It's easy to show that 

= -VsOr)F(t), = - nj H s (x)5(t). (16) 

Hereinafter we will consider the generalized functions 

u(x, t) = u(x, t)Hp(x, t), G(x, t) = G(x, t)H s (x) H(t), 

where u(x,t) is the classical solution of BVP. 

Let us define the action of wave operator on u(x,t) . Since [u]d = —u, executing generalized 
differentiation and using Eq. (|TJ) in the field of differentiability, we get 



du 

c~ 2 Hg (x) uq (x) 5(t) — c~ 2 Hg (x) mq (x) S(t) + G, 



D c u (x, t) = -—5 S (x) H{t)-H (t) {un 3 5 s (x)) , 3 



where (3 (x,t) 5s(x)H(t) is the simple layer on {S x R + } , S(t) is Dirak function , |^ = u^rii 
is a derivative along normal vector n to S . Notice that density of simple and double layers are 
here defined by use of boundary conditions (a part from which, in depend on solved BVP, are 
known) and by given initial conditions. 

We can present the solution of the Eq. (TTTj) as the convolution of its right part with the 
fundamental solution U (x, t) : 

U C U = 5(x)5(t), (18) 

U = for t < and ||x|| > ct (radiation conditions) . (19) 
We shall name it Green function of Eq. (JTJ. 
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The solution (fT6|) is represented in the manner of the following convolution (*): 

°" V(.r) /i it) - [i * >ni t d s (.r) H (/) ] ., - 

(20) 



u ( x , t)=u (x, t) Hg (x) H (t) = -U (x, t) * -^-5 S (x) H (t) - (u * un^s (x) H (tfj 
-cT 2 (u * H s (x) u (x)) , t -c- 2 U * H s (x) u (x) + G * U, 

\ X J X 



where symbol " * "means that convolution on x is taken only over x . Moreover the solution 

X 

(poll is single solution in the class of functions, allowing the convolution with U (x, t) . 



If initial data and G are null then 

du 
On 



u (x, t) = -U (x, t) * —5 S (x) H (t) - U,j *urij5s (x) H (t) . (21) 



This formula expresses the BVP solution through boundary values of function u and its 
normal derivative and similar to Green formula for solutions of Laplace equation [2]. Howev- 
er, on the strength of particularities of the fundamental solutions of hyperbolic equations on 
wave fronts, which type depends on space dimensionality , its integral representation gives the 
dispersing integrals (in second summand). 

For building its regular integral presentation we introduce the antiderivatives on t functions: 

W (x, t) = U (x, t)*8 (x) H(t) = U (x, t)*H (t) d t W (x, t)=U (x, t) ; (22) 

N dW(x,t) dW(x,t) 
H(x,m,t) = — - 



1 " v 1 " I r\ • • -I r% 

oxi am 

dm 



Easy to see that they also are solutions of Eq.(JH) when G = H(t)5(x) and G = H(t) 9S g^ 
accordingly. 

Theorem3.1. In D'(R N+1 ) BVP solution satisfies to the equation: 

u (x, t) = —U (x, t) * ^5 S (x) H (t) - (w n *u nj 5s (x) H (t)) - (23) 
— W,j * uo (x) rij (x) 5s (x) — c~ 2 U * HZ (x) u (x) — cT 2 ( U * H s (x) u (x) ) , t +G * Uyu 

X X \ X J 

Proof. It's easy to show, using determination of derivative of generalized functions and 
continuity u , that 

(urijSs (x) H (t)) , t = u(x, t)rij(x)5s (x) H (t) + u(x, 0)rij(x)8s (x) 5 (t) . 
Using this relation, last formula (|2*2T) and properties of convolution differentiation , we have 

U * urij5 s (x) H (t) ) ,j = W,j *u(x, t)rij(x)5s (x) H (t) + W,j * Uo{x)nj(x)5 s (x) 



Substituting these correlations in formula (|27Hl . we get the formula of the theorem. 

From theorem follow that solution of the problem is completely defined by boundary value 
of the normal derivatives of u (x, t) and its velocities u . In analogy with presentation of the 
solutions of the Laplace equation, these formulas possible to name dynamic analogue of Green 
formula. 

The formula (|2Tfl) possesses the advantage in comparison with ([271 . since it allows immediately 
to go to its integral record without regularization of integral function on front, as it was earlier 
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offered in paper [3]. For x G S formula (f23|) gives, what we show hereinafter, boundary integral 
equation for BVP solution. If one of the boundary function is known, then after solving BIE 
on border, we can find second boundary function. Whereupon formula (f23j) gives the solution 
u(x, t) in D~ . 

4. Dynamic analogue of Gauss formula. We introduce the functions 

U (x, y,t) = U(x-y,t), W (x, y,t) = W (x - y,t) , H (x, y, m,t) = H (x - y, m, t) , 

which, on the strength of properties of symmetries of the wave operator and 5 -functions, satisfy 
to following symmetries correlations: 

dW dW 

U{x,y,t) = U{y,x,t), W{x,y,t) = W{y,x,t), -— = --—, 

oxj dyj (24) 

H (x, y, m, t) = -H (y, x, m, t) = -H (x, y, -m, t) . 
Lemma 4.1. In D'(R N+1 ) the dynamic analogue of Gauss formula has the form: 

-W,i * Hi (x) 5 (x) - c- 2 (U * H s (x)) , t = R- D (x, t) (25) 

X \ X / 

P r o o f. If in both parts of Eq. (JTK| for U the convolution with Hp(x, t) are taken, after using 
a property of convolution differentiation and formula (jTfil) we get: 

-U, d *n,(x)5s (x) H(t) - c~ 2 (U, t *Hg (x) 6{t)) = H D (x, t) 

With regard for t|22j) , flipping differentiation on t in the first summand and executing convolution 
on t in second one, we get the formula of this lemma. Integral record of this formula depends 
on dimensionality of the equation. 

The formula (|25|) is an analogue of the known formula of Gauss for double layer potential 
(p. 406, [2]), which gives the integral record of characteristic function of a set with use the 
fundamental solution of Laplace equation. The Gauss formula is often used for building the 
boundary integral equation of boundary value problems for elliptical equations and systems. 
Similarly it's possible to use the dynamic analogue of Gauss formula for constructing BIE in 
hyperbolic case. But here the dynamic analogue of Green formula is used for building BIE of 
setting problems. 

Hereinafter we shall give integral representation of formulas of the theorem 3.1. and lemma 
4.1 for space of dimensionality N=l,2,3, that is most typical for problems of mathematical 
physics. 

5. BIE for plane BVP. When iV = 2 we have a plane problem. At first let's consider the 
problem with null initial conditions. 

We denote as dS(y) the arc length differential on S in point y , S t (x) = {y G S : r < ct} , 
s t( x ) = {y £ S~ : r < ct}, r=\\x- y\\ , dV(y) = dy x dy 2 ■ 

Theorem 5.1. By N = 2 function u (x, t) with null initial conditions (uq = 0, Uq = 0) , 
has following integral representation: 

2tt J J V Mv) r My) J Vc 2 r 2 - r 2 v 1 



Sr(x) 
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Or, if to change the integration order, it has the form 



u{x,t) = J dS [y) J 



t 

du (y, t — t) dr 



St(x) r/c 



dn (y) ^( c 2 r 2 - r 2 ) 



+ c_ I \dr dS{y) fru(y,t-r)dr 



2tt J rdn(y) J ^/( c 2 t 2 - r 2 ) ' 

St(x) r/c 

For x G S second integral on the right part is singular and taken in value principle sense. 
Proof. For N = 2 Green function has a following type ([2], p. 206) 



Calculating from formula (|22|) we determine 



^>2^— 2 /y 2 



U = -^M=*r> R=J.ri + 4- (27) 



rtr/ \ H(ct-R), f ct- VcH 2 - R 2 \ . ctH(ct-R)x j m j . , 

W(x,t) = K — ^ln V — , H{x,m,t) = - - -. \„ ' ^T- (28) 



2tt V R J 2tt VcH 2 - R 2 R 2 



If to write the convolutions (|23|) in integral type with regard for these relations, that we 
get the formula of the theorem. Notice that for x S integrals on the right are regular and 
therefore for such x on the right and on the left the regular functions stand . We shall prove 
formula ((SB)) also for x G S with regard for determinations Hg(x) . 

We denote e -vicinity of the point x (e <C ct,t > 0) through III e (x) = {y : r < e} , 
S~ (x) = S -U1 £ (x) , S+ (x) = S + -U1 £ (x) , o £ {x) = {y G S : r < e} , S £ = S - £ , in; = 
S- n ni £ , Ult = S + n UI £ , Vf (x) = {yeS ± :r = e}, r, } = |-. 

Let x* G S . We transform the sidebar S in vicinity of the point x* , avoiding it on e - 
semicircle in S~ (e << ct,t > 0) . From dynamic analogue of the Green formula for sidebar 
S £ + T~ we have in point x = x* 

t 

c f v , al . f du{y,t-r) dr 



= ^ J H{ct-r)dS (y) J 

S £ {x*)+r-(x*) r/c 



dn (y) y/ 



^>2^— 2 ty 2 



+^ / H(ct-r) 1 -^-dS(y) tfTil{y ' t - T)dT 



2vr J r dn(y) J \f 



c 2 r 2 — r 2 



S^x*) r/c 



t 



T~{x*) r/c 



dr f ru(y,t-T)dr 



c 2 t 2 — r 2 



For e — > +0 first integral, on the strength of weak singularity of the integral functions, strives to 
integral on S , second one tends to integral in the value principle sense, which also exists, since 



integral function has a singularity of the type r 1 and contains function = ^(y)^—^- 

dr _ „ f„\ (Vj-Xj) 



dn(y) 

which for y — > x asymptotically equivalent to = rij (x) yyj , antisymmetric relative to 

point x . 
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Let consider the last summand, denoted J £ (x) . On T~ there are r = e, = —1, dS (y) = 
e d9 , where 6 is an polar corner with top in point x ; Q\ and 62 are the corners of endpoints 
T~ , numbered in order at pass-by of the sidebar T~ in positive direction. With regard for these 
relations 

Tu(y,t — T)dr (62 — Q\) c f ru(y,t — r) dr 



j F (x) = — / -de 



2 /y-i 2 



2tt J e J y/c 2 r 2 — r 2 27r J yj c *t* — r' 

81 e/c e/c 

t 

lim (# 2 — #1) = — 7r, lim J £ (x) = —0, 5 / u (y,t — r) dr = —0, 5m (x, t). 

e^O e— >0 J 



Passing this summand in the left part, with regard for determinations: H^(x,t) = 0,5 for 
x G S , we get the formula of theorem on the border also. Since on the left and on the right in 
Eq. |26| the regular generalized functions stand, on the strength of Dyubua-Reymon lemma ([2], 
p. 97) this equality, which equitable in generalized functions class , also equitable in classical 
sense. This theorem has been proved. 

By solving the BVP the Eq. (|26|) gives the boundary integral equation for determination of 
unknown boundary function: 

nu s (x,t) = 

t t 
du(y,t — r) dr f 1 dr f TUs(y,t — T)dT 



V J J dn{y) ^Jr 2 - r 2 /c 2 J rdn{y) yy, J ^r 2 - r 2 /e 

St(x) r/c St(x) r/c 

In the case of the first BVP the left part of this equation and the second integral on the right 
are known, they are defined with use of boundary condition, but the first integral contains the 
kernel with weak singularity on front of the Green function. Solving it we define the normal 
derivative of u on boundary whereupon formula (|26jl allows to compute the solution in any 
point. 

In the case of the second BVP we have singular BIE for determination of unknown boundary 
values of u . After its solving we can define u(x, t) on boundary whereupon formula l|2fij) defines 
the solution fully. 

In the case of nonzero initial conditions the following theorem gives the solution of problems. 
Theorem 5.2. By N = 2 BVP solution has the following integral representation: 

9 - ) a f f du(y,t-r) 1 dr \ dS(y) , 

+ d_ r u Q (y)dV(y) | r u Q (y)dV{y) { c f dr f G{y,t-r)dV(y) 



dt J cV cH 2 — r 2 J cV c 2 t 2 — r 2 J J V' 

St(x) S t -(s) S~(x) 

ct dr dS (y) 



1 ^-t2i '^■'21 ^2 



r dn (y) y/ c 2 r 2 - r 2 



St(x) 



Proof follows from theorem 4.1. and 5.2. if to write the convolution with initial data in integral 
type. Here the integrals from the second before the fourth ones comply with Poisson formula for 
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Caushy problem . The last additional summand with initial data is conditioned presence of the 
boundary. 

Using relations (|2"TJ) . dynamic analogue of Gauss formula possible to be write in integral type. 
L e m m a 5.1. By N = 2 dynamic analogue of Gauss formula has the following type: 

Sr(x) V VI ' S-(X) 

where integral in the value principle sense is taken for boundary points . 

Proof. Formula of the lemma 4.1 , with regard for (j2ZJ) , J2H| , possible to write so 

J VcH 2 -r 2 rdn (y) Kyj cdt J ^JcH 2 - r 2 S\ ) \ ) \ ) 

St(x) St(x) 

From this formula by elementary transformations we get the formula of the lemma. 

We shall show that this equality, equitable in the field of regularity , is saved also for x G S 
if the singular integral on the left, which contains the strong singularity on r , is taken in the 
value principle sense. 

Similarly (|29jl for area without and with e -vicinity of the point x we get 



ctH(ct-r)l dr . . d f H (ct — r) dV (y) 



Vc 2 r 2 - r 2 rdn(y) cdt J y ' c 2 t 2 - r 2 

ctH(ct — r)l dr . . d f H (ct — r) dv (y) 

— --- — r-rdb (y) H — — / : — = 2%. 

Vc 2 r 2 - r 2 r dn (y) vy; cdt J Vc 2 r 2 - r 2 

5 e +r+ sr(x)+w e 



(30) 



Under e < ct integrals on T e , III £ are easy calculated by transition to polar coordinate system. 
On = ±1 consequently to the sign. Mark that 

2tt e 

H {ct — r) dV (y) f f rdr 



I £ (x,t)= / \ / = dff — — = 2vr(ct - VcV 2 - r 2 



\/c 2 t 2 — r 2 J J yc 2 r 2 - r : 

m e (x) o o 

dl ( ct 
lim — - = 2ttc lim 1 , 1 = for Vt > 0. 



e- 



+ dt e^O ^ y/ t 



c 2 t 2 — r' 



If both equalities in (f3n|l to add and divide on 2, with regard for properties of the symmetries of 
integral function (f2~i|) and going to limit on e — > , then 

-V.P. / —===-—— dS(y ) + — / - Ky) = tt. 

7 Vc 2 r 2 -r 2 rdn (y) cdt J yjc 2 r 2 - r 2 

St(x) St{x) 



So, with regard for determinations H s (x) , formula (|29|) is equitable for any x 
6. BIE of BVP for N=3. 
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Theorem 6.1. For N = 3 BVP solution is represented in the form: 

Anu(xt)= I f l du (y> t - r / c ) | u(y^- r / c ) g ln r \ dS ( ) 

J \r dn{y) c dn(y) J 

St(x) 

-V.P. J u { y, t -r/c)^l iS{ y) + ^l J M y)^S(y) + 

St(x) St(x) 



r=ct r=ct 



Proof. For AT = 3 Green function U (x, t) is the double layer on the cone K t = 
{(x, t) : \\x\\ = ct} : 

frfo*) = - * ( * 4 "^ /c) , ^=v^+^+^i- ( 31 ) 

For any (x, t) & D (i? 4 ) it defines the linear functional: 
After calculation on formula (l22l) we find 



For building the integral representation of dynamic analogue of Green and Gauss formulas 
we use the following equalities, which be get using the determination of generalized function 
convolution : 



a(x)5(t-R/c)*f(x,t)H(t) = H(t) j a (x - y) f (y,t - r/c) dS (y) 

St(x) 
t 

P(x,t)5 s (x)H(t)*f(x,t)H(t) = H(t) J dr J (3(y,t-r)f(x-y,T)dS(y) 

St Or) 



a(x)5(t- R/c)* f3(x,t)5 s (x)H(t) = H(t) J a(x - y) p (y,t- rj c) dS (y), 

St(x) 

d f 

a (x) 8(t- R/c) * 7 (x) 5 s (x) = — / a (x - y) 7 (y) dS (y) 

St{x) 

a{x)5{t-R/c)*H s {x)=c- 1 H{t) J a {x - y) Hg {y) dS (y). 

r=ct 
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We compute the convolution in Eg. ()23|) with regard for these correlations. First summand in 
this equation is 

r \ i x du . , , „ , s H (t) f 1 du (y, t — r/c) , _ . . 

- u{x , t )*- Ss ( X )H( t) = -n J - 

St(«) 

The second summand ( — 4tt W,j *u(x, t)nj(x)5s (x) H (t) ) is an amount of two convolutions: 
' Xj S(t- \\x\\ /c) *u(x,t)n j (x)5 s (x)H (t) = c~^H(t) [ u(y,t-r/c)^^rdS(y), 



St(x) 

-%H(ct- \\x\\)*u{x } t)n j {x)S s (x)H(t) = 



y ""^^^-^ y «(v,f)^ds(y) 

5t(x) " St(x) 



d f dlnr f dr 



,-i 



Third summand 

-4tt * u (x) nj (x) 5 S ( x ) = ^ J ^yj u ° (f ) rf,S (f) _ y M ° (y) rfS 

St(x) " S t (x) 

The last three summand give the known Kirhoff formula for solution of the Caushy problem [ 
1,2] with initial conditions. Summing these convolutions, we get the formula of the theorem. 

In the formula of the theorem in the part, depending on initial data, there is one summand 
conditioned by presence of the boundary S . It disappears for t > t*(x), t*(x) = max ^ x ~ y ^ since 

yES 

S t (x) = S and integral does not depend on t . 

In formula (fT6|) for x S, t > all integrals exist. Its proof for x G S like to plane case. 
Herewith the strong singularity has the second summand on the right. In this case 

hm J u(y,t- ~) \ dS (y) = _ J^~2 J u (y^~ ^) dS (y) = ~ 2nu > 

r- r- 

It's clear that for x G S formula (|3Tf saves kind if corresponding singular integral to calculate 
in the value principle sense and take into account the value Hg (x) on S . 

For x G S this formula gives the BIE for problems solving , it is singular for the second 
BVP. 

Using correlations (fT9|) - ([22|) the dynamic analogue of Gauss formula also possible to write in 
integral form. 

L e m m a 6.1. By N=3 

St(x) [r=ct St(x) ) 

When t > t* (x) thence the known Gauss formula [ 2] follows: 

/sfcj® <»> = (33) 

s 
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Proof of this formulas similarly is like to proof of the lemma 5.1. Notice that for N = 2 
two-dimensional analogue of the Gauss formula is not followed from lemmas 5.1. 

7. Solution of the BVP for Dalamber equation (N=l). 

Theorem 7.1. For N = 1 BVP solution has the following integral presentation 

t t 

2u = cH {ct — \x — 02]) J u, x (a 2 ,r) dr — cH {ct — \x — a x \) J u, x (a x , r) dr+ 



\x—a,2\/c \x— a\\/c 

\x — a x \ 



+ sgn (x — a x ) H (ct — \x — a x \) u I a x ,t 



c 



— sgn (x — a 2 ) H (ct — \x — a 2 \) u \a 2 ,t J + 

+c _1 J ii (y) H(ct — \x — y\)dy + u (x + ct) H$ (x + ct) + u (x — ct) Hg (x — ct) . 

ai 

Proof. Denote xx — x . In this case ([2], p. 206) 

U (x, t) = -°-H (ct-\x\), U, t = - C -b (t-\x\ /c) , 

W (x, t) = ~H (ct - \x\) (ct -\x\), W, x (x, t) = (ct - \x\) sgn x, 



(34) 



(35) 



1, x > 0; 

sgn a; = { 0, x = 0; (36) 
-1, x<0. 

In this case it is impossible to use the formula (|23|) for building of the integral analogue of the 
Green formula , because in it some functions are not determined . 

It is possible to get the similar formula, if put u as zero outside of given interval and 
considering an action of wave operator on it in the class of generalized functions. We enter 
otherwise to use the formula ([20"|) . We increase the domain of definition u(x,t) in the band 
in R 2 x R + : {a x ^ X\ ^ a 2 , —00 < x 2 < 00, t > 0} . Then boundary S will consist of 
two direct lines X\ — Oi, X\ — 0,2 ) 

which external normals have coordinates (-1,0) and (1,0) 

accordingly, = raij^- for x G S , Hg (x) = H (x x - a x ) H (a 2 - x x ) , n x 5 s (x) = = 
—5 (x x — a x ) + 5 (x x — a 2 ) . The formula (|23| of the theorem 3.1 is converted to type: 



u = U 2 * -r^H (t) (8 (x -a 2 )-5(x- a x )) + * u (x, t) H (t) (5 (x - a 2 ) - 5 (x - a x )) + 
ox ox 

+c~ 2 (u 2 , t * u (x) H s (x) \ + c~ 2 U 2 * G. 

(37) 

Here, all convolutions are taken with Green function for N = 2 and its antiderivative. On base 
of the method of the lowering on x 2 , rolling up on x 2 , since u does not depend on x 2 , we get 

n du (a 2 , t) rr , , - . . du (a x ,t) , . 

H = U(x- a 2 , t) * (t) -U(x-a x ,t)* \r^H (t) + 

t ox t ox 

+W, X (x - a 2 , t)*u (a 2 , t) H (t) - W, x (x-a x ,t)*u (a x , t) H (t) + W, x (x - a 2 , t) u (a 2 ) - 
—W, x (x — a x , t) u (a x ) — c~ 2 ( U * u (x) He (x) ) — c~ 2 U, t * u (x) HZ~ (x) + U * G. 

\ X J X 

(38) 
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Substituting (|36|) into (j38|) and executing integration, we get formula (J3IJ) . 

Easy to show that it is equitable also for x = ai, x = a 2 (with regard for pfij) ). For this 
it is enough to write formula (jUj) for interval (a'x,a 2 ) = (ax +£,02) ((ax,a2 — e) ). Supposing 
x = a x and e — > +0 we have: 

If f 1 

= lime < (ct — d) / u, x (a 2 ,r) dr — H (ct — e) / it )a . (ai, r) dr > + 



-1 



+i? (ct — d)u (a 2 , t — d/c) + sgn (— e) (ct) it (ai, t) + c y -it (?/) i?(ct — |ai — y|)djH- 

ai+e 

+m (ai + ct) Hg (a x + ct) + w (a x — ct) i?^ (ai — ct) = 

r t « 1 

= c < H (t — d/c) J u, x (a 2 , t) dr — J u, x (ai, r) dr > + H (t — d/c) u ^a 2 , t ^ + 



c 

a> 



—H (t) u (ai,t) + c 1 J iio(y) H(t — |a x — y\ /c)dy + u (a x + ct) (a x + ct) + 

ai 

+u (ai — ct) (a x — ct) 

(d = |ai — a 2 |). Transferring the summand —H(t)u(ai,t) into the left part, with regard 
for values of the characteristic function on border, we get the formula of the theorem for left 
endpoint. By the similar way the formula is proved for x = a 2 . As a result on the end of 
the interval (ax, a 2 ) we have the following equations for determination of unknown boundary 
functions: 

t t 
u(a±, t) = cH (ct — d) J u :X (a 2 , r) dr — cH(t) J u iX (a±, r) dr + H (ct — d)u (a 2 , t J + 

d/c 

(12 

+c _1 J u (y) H(ct — |ai — y\)dy + u (a\ + ct) H (d — ct) H (t) /yra x = a x ; 



O-l 

t 



u(a 2 t) = cH (t) J u, x (a 2 , t) dr — cH (ct — d) j u, x (a 1; r) dr + H (ct — d) u ^a 1; t j + 



d/c 

02 

-1 



+c / u (y) H(ct — \a 2 — y\)dy + u (a 2 — ct) H (d — ct) H (ct) pjin x = a 2 . 



Under given u, x (a^, t) , k = 1, 2 , we get two functional equations with lagging argument for 
determination u on boundary of the area, which can be solved incremental on time from t = . 
Under the known u(a\,t) , u (a 2 ,t) we have the system of two integral equations. 
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Conclusion. Using this method it is possible to build the similar formulas and boundary 
integral equations for solutions of BVP in space of greater dimension ( N > 3 ). It is particularly 
efficient in BVP for systems of the equations of mathematical physics, when the Green matrix 
of system can me constructed. Herewith the type of the equations is unessential , it can be also 
elliptical, as, for instance in problem of stationary diffraction of electromagnetic waves [5], or in 
problem of the elasticity theory [6], parabolic or mixed type in problem of thermoelastodynamics 
[7]. But particularly efficient this method for solving of the hyperbolic equations, where using 
classical methods more difficult, but sometimes and simply impossible (refer to, for instance, 
[6,8,9]). 

The study of solubility of built BIEs presents the independent problems of the functional 
analysis since these equations do not pertain to well studied classical ones. However we notice 
that using the computing methods on base of the methods of boundary element with transition 
to discrete analogue BIE, allows effectively to build the solutions of like problems [10]. 
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Alexeyeva L.A. Nonstationary boundary value problems for wave equation and their 
generalized solutions 

Abstract. The multivariate analogue of Dalamber's equation in the space of generalized functions is 
considered. The method of generalized functions for the building of solutions of nonstationary boundary 
value problems for wave equations in spaces of different dimensions is elaborated. Dynamic analogues 
of Green and Gauss formulas for solutions of wave equation in the space of generalized functions are 
built. Their regular integral representations and singular boundary integral equations for solving the 
nonstationary problems are constructed for the spaces of the dimensions 1,2,3. The method of obtaining 
of conditions on fronts of shock waves is stated. 
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